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ABSTRACT 


We have discovered kiloparsec-scale extended radio emission in three narrow-line Seyfert 1 galaxies (NLS1 s) 
in sub-arcsecond resolution 9 GHz images from the Karl G. Jansky Very Large Array (VLA). We find all 
sources show two-sided, mildly core-dominated jet structures with diffuse lobes dominated by termination 
hotspots. These span 20-70 kpc with morphologies reminiscent of ER II radio galaxies, while the extended 
radio luminosities are intermediate between ER I and ER II sources. In two cases the structure is linear, while 
a 45° bend is apparent in the third. Very Long Baseline Array images at 7.6 GHz reveal parsec-scale jet 
structures, in two cases with extended structure aligned with the inner regions of the kiloparsec-scale jets. 
Based on this alignment, the ratio of the radio core luminosity to the optical luminosity, the jet/counter-jet 
intensity and extension length ratios, and moderate core brightness temperatures (< 10*** K), we conclude these 
jets are mildly relativistic (J3 < 0.3, 6 ~ 1-1.5) and aligned at moderately small angles to the line of sight 
(10-15°). The derived kinematic ages of ~ 10^-10^ y are much younger than radio galaxies but comparable 
to other NLS Is. Our results increase the number of radio-loud NLS Is with known kiloparsec-scale extensions 
from seven to ten and suggest that such extended emission may be common, at least among the brightest of 
these sources. 

Subject headings: galaxies: active — galaxies: jets — galaxies: Seyfert — galaxies: individual (J0953H-2836, 
J1435H-3131, J1722-I-5654) — radio continuum: galaxies 


1. INTRODUCTION 

Narrow-line Seyfert 1 galaxies (NLSIs) are a peculiar 
class of active galactic nuclei (AGN). They resemble ordi¬ 
nary Seyfert 1 galaxies but their optical broad lines are nar¬ 
rower (EWHM(//p) < 2000 kms“'). The black hole masses 
of NLS Is are relatively small (10^ - 10* M©; e.g., Zhou et al. 
2006) and their accretion rates are high, often a substantial 
fraction of the Eddington rate (Boiler et al. 1996). It has, how¬ 
ever, been suggested that this is a result of systematic under¬ 
estimation of the black hole masses (Decarli et al. 2008; Mar¬ 
coni et al. 2008; Calderone et al. 2013). About 7% of known 
NLS Is are radio-loud (Komossa et al. 2006) and it is clear 
that jets are present in some of these (e.g., Zhou et al. 2003; 
Yuan et al. 2008). The detection of gamma-ray emission from 
several radio-loud NLS 1 s provides evidence that these jets are 
relativistic (e.g., Abdo et al. 2009). This was unexpected be¬ 
cause jets in AGN are typically associated with large black 
holes, while high accretion rates in stellar-mass black hole bi¬ 
naries normally correspond to quenched-jet states (e.g.. Boro- 
son 2002; Maccarone et al. 2003). Study of the jets in NLS Is 
will therefore help us understand how the processes underly¬ 
ing jet production scale with black hole mass (e.g., Heinz & 
Sunyaev 2003; Eoschini 2012). 

Radio emission on kiloparsec scales has been found in only 
a few NLS Is. Together, Whalen et al. (2006), Anton et al. 
(2008), Gliozzi et al. (2010), and Doi et al. (2012) reported six 
examples from the 1.4 GHz Eaint Images of the Radio Sky at 
Twenty centimeters (EIRST) survey (Becker et al. 1995) and 
4.8 GHz Australian Telescope Compact Array (ATCA) ob¬ 
servations. The resolution of these observations was limited, 
however. The 5" EIRST beam only resolves sources >2", 
while the ATCA beam was about a factor of two better. NLS 1 s 
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are known out to z = 0.8 so this corresponds to >10 kpc, leav¬ 
ing a range of kiloparsec scales largely unexplored. Recently, 
Doi et al. (2014) reported extended emission in a radio-quiet 
NLSl. 

In this work, we present the detection of kiloparsec- 
scale radio emission in three radio-loud NLSls (J0953H-2836, 
J1435H-3131, and J1722-I-5654). These are the first such de¬ 
tections in these sources, all of which were detected but un¬ 
resolved by the EIRST survey (White et al. 1997). This in¬ 
creases the number of radio-loud NLSls known to be ex¬ 
tended on kiloparsec scales from seven to ten. We also present 
parsec-scale radio images showing bright jet-like structures. 
Where necessary, we have assumed a flat A-Cold Dark Matter 
cosmology with Da = 0.7 and Hq - 70 kms“*Mpc“^ Spec¬ 
tral indices are specified with flux density 5y oc v“ at fre¬ 
quency V. Velocities are quoted as fractions of the speed of 
light, p - vjc. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed three radio-loud NLSls with the Karl G. Jan¬ 
sky Very Large Array (VLA) and the Very Long Baseline Ar¬ 
ray (VLBA). The VLA observations were made to obtain ra¬ 
dio coordinates accurate enough for correlation as part of a 
VLBA program monitoring the 15 known radio-loud NLSls 
above 0° declination with archival centimeter-band flux densi¬ 
ties exceeding 30 mJy. (Richards et al., in prep.). Among the 
15, only these three lacked sufficiently precise coordinates. 
Properties of our targets and our refined coordinates are listed 
in Table 1. 

The VLA observations were made on 15 Leb 2014 and con¬ 
sisted of a single 10 minute snapshot of each source sand¬ 
wiched between short observations of a nearby phase calibra¬ 
tor with a known position (better than 2 mas). A short scan on 
3C 286 was used to calibrate the flux density scale, applying 
the appropriate model to account for resolved structure. At the 
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time of observation, the VLA was transitioning from BnA to 
A configuration, with 27 active antennas on baselines ranging 
from 0.3 to 31.5 km. The X-band receiver provided 2 GHz 
bandwidth centered on 9 GHz, which was sampled with the 
8-bit sampler with 2 s integration time. Requantizer gains 
were applied and the data were converted to UVFITS format 
using CASA (McMullin et al. 2007). The data were then cali¬ 
brated in AIPS following the standard procedures in the AIRS 
Cookbook and deconvolved, self-calibrated, and imaged in 
DIFMAP (Shepherd et al. 1994). After initial calibration of 
the phase calibrators, cycles of phase-only self-calibration and 
CLEAN were repeated until a few times the residual root- 
mean-square (rms) noise level was reached. Amplitude self¬ 
calibration was then applied to determine the final complex 
gains. To measure source positions, the gains derived from 
the phase calibrators were interpolated onto the science obser¬ 
vations. A circular Gaussian component model of the core re¬ 
gion was fitted to the visibilities and the position of the bright¬ 
est component is reported in Table 1, accurate to about 0.1 
beam, or about 15 mas. For imaging, the science target data 
were cleaned and self-calibrated using the same procedures 
used for the phase calibrators. 

VLBA observations were made on 8 Feb 2014 using the 
wideband C-band receiver configured to provide 256 MHz 
total bandwidth in 8 sub-bands equally spaced over the 
480 MHz band centered on 7.632 GHz with good weather 
at all 10 antennas. The total integration time, listed in Ta¬ 
ble 1, was split into scans separated in time to provide good 
(u,v) coverage. In each scan, 49 s of full polarization data 
were recorded as 2-bit samples for a total 2 Gbps data rate, 
with science target observations interleaved with observations 
of a nearby phase reference. Observations at 4.9, 15.4, and 
23.8 GHz were also carried out, but we report only 7.632 GHz 
total intensity results here. Full results from the VLBA pro¬ 
gram, including polarization measurements and the other fre¬ 
quencies, will be reported elsewhere. After correlation with 
the DiFX correlator (Deller et al. 2007) using our refined coor¬ 
dinates, the phase-referenced VLBA data were processed fol¬ 
lowing the AIPS Cookbook, applying gain solutions from the 
phase references to the science targets. A CLEAN model was 
constructed, then phase-only self-calibration with a 10 minute 
solution interval was applied and a final CLEAN was per¬ 
formed. We verified our flux density scale, accurate to about 
5%, by comparing several calibrators with preliminary results 
from Effelsberg 100 m telescope observations made by the 
Fermi-GST AGN Multi-frequency Monitoring Alliance (E- 
GAMMA; Euhrmann et al. 2007; Angelakis et al. 2010). 

3. RESULTS 

The maps we obtained from our observations are shown 
in Eigures 1, 2, and 3 and Table 2 lists quantities computed 
from our results. Using 1.4 GHz EIRST flux densities and 
our 9 GHz VLA flux densities, we have determined the non- 
simultaneous spectral index a. The formal uncertainty is 
+0.04, assuming 5% uncertainties for the flux densities, not 
accounting for non-simultaneity. These spectral indices are 
lower limits because our higher resolution maps may resolve 
out more extended emission than in the EIRST data. 

On kiloparsec scales, all of these sources show similar mor¬ 
phologies. The emission is dominated by a bright core strad¬ 
dled by a pair of extended lobes. These lobes are edge- 
brightened, terminating in hotspots. Like in ERII radio galax¬ 
ies, the ratios of the hotspot separations to the full extents of 
the sources exceed 0.5 (Eanaroff & Riley 1974). Extended 


emission between the bright core and the brighter of the two 
lobes is found in all sources, while the fainter lobe is rela¬ 
tively isolated. On parsec scales, we find an isolated core in 
J1435+3131 while the others show extension roughly aligned 
with the kiloparsec-scale structure. 

3.1. J0953+2836 

The emission in the tapered VLA map of J0953+2836 (Eig- 
ure 1) spans 1071 (measured at the second contour, follow¬ 
ing Scheuer (1995)), a projected length of 70.2 kpc. The 
brighter lobe lies along position angle (PA) 212°, measured 
north through east. Prom the VLBA map, the parsec-scale 
image shows extended emission with a ridge along PA 232°. 
The ridge brightens slightly toward the SW. A two-component 
Gaussian model fits the visibilities well, with the SW compo¬ 
nent brighter by a factor of 1.2. Identification of the parsec- 
scale core is uncertain and it is unclear whether this is a one- 
or two-sided jet. 

3.2. J1435+3131 

The kiloparsec-scale emission in J1435+3131 (Pigure 2) 
spans 5700, or a projected length of 30.6 kpc. A sharp bend 
is evident between the core and the NE kiloparsec-scale lobe. 
A narrow, continuous trail of emission exits northward from 
the core, then bends about 45° to the northeast, terminating in 
a faint knot. A straight line extended from this knot along the 
direction after the bend is closely aligned with the brighter of 
two hotspots in the disconnected radio lobe. The parsec-scale 
image shows only a compact core. No kiloparsec-scale emis¬ 
sion is evident between the core and the fainter S W radio lobe, 
which has a single hotspot. 

3.3. J1722+5654 

On kiloparsec scales, J1722+5654 (Pigure 3) contains a 
bright central core straddled by a linear structure along PA 
47°. The emission spans 4708 which corresponds to a 
22.8 kpc projected length. The parsec-scale structure is also 
linear, with three apparent emission regions along PA 60°. 
The brightest parsec-scale emission is from the SW compo¬ 
nent; if this is identified as the core, then the parsec-scale 
morphology is a one-sided jet directed toward the brighter 
kiloparsec-scale lobe. This component is 2.3 times brighter 
than the NE component. We note that the structure could also 
be interpreted as a two-sided jet with a faint core at the central 
emission region, in which case the brighter ejected component 
is directed toward the fainter kiloparsec-scale lobe. Misalign¬ 
ments as large as 180° between the parsec- and kiloparsec- 
scale jets have been seen (Kharb et al. 2010). However, the 
brightness temperature of the SW component is more than 
20 times higher than either other component, which favors its 
identification as the core (see Section 4.3). 

4. DISCUSSION 
4.1. Morphologies 

The parsec- and kiloparsec-scale morphologies of these 
three sources are similar to those found in other jetted AGN. 
Their kiloparsec scale morphologies resemble those found in 
blazars, radio galaxies, and other radio-loud NLSls (Cooper 
et al. 2007; Kharb et al. 2010; Doi et al. 2012). Extended 
two-sided emission is present in all three cases, with one lobe 
substantially brighter than the other, presumably because it is 
approaching us. These resemble ER II radio galaxies, with 
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TABLE 1 

Source properties and refined J2000 coordinates. 


Source 

NVSS Name 

Z 

Mbh 

Sl.4 

Mv 

R.A. 

Decl. 

Phase Reference 

VLBA Time 




(lO^Mo) 

(mJy) 


(h m s) 

(O , 

VLA / VLBA 

(min) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

J0953+2836 

J095317+283601 

0.659 

6.3 

47.94 

-24.04 

09 53 17.102 

28 36 01.559 

J0956+2515/J0954+2639 

4.8 

J1435-I-3131 

J143509+313149 

0.502 

3.2 

44.72 

-23.01 

14 35 09.495 

31 31 47.864 

J1435+3012/J1422+3223 

9.8 

J1722+5654 

J172206+565452 

0.426 

2.5 

39.83 

-23.58 

17 22 06.029 

56 54 51.696 

J1657+5705/J1727+5510 

13.0 


Note. — Columns are as follows: (1) lAU name (J2000); (2) corresponding source in the NRAO VLA Sky Survey (NVSS, Condon et al. 1998); (3) redshift 
from Hewett & Wild (2010); (4) black hole mass from Yuan et al. (2008); (5) integrated 1.4 GHz flux density from White et al. (1997); (6) absolute V-band 
magnitude computed from Sloan Digital Sky Survey Data Release 10 (Ahn et al. 2014) g and r magnitudes using V = ^-0.52(^-r)-0.03, after Jester et al. 
(2005) with no correction for galactic extinction; (7)-(8) new J2000 coordinates determined from our VLA observations with an uncertainty of ± 15 mas; (9) lAU 
name (J2000) of phase calibrators used for VLA and VLBA observations; (10) total VLBA integration time. 


TABLE 2 

Measured properties. 


Source 

^proj 

f^kpc 

core,9 

„kpc 
tot,9 

7 -hs,j 

9 

a 

rext 

^1.4 

Rb. 

logRv 

(/3cose)hs 

(/ScosS)^ 

Age 

core,7.6 

cPC 

^ tot,7.6 

Tb 


(kpc) 

(mJy) 

(mJy) 

(p,Jy beam“^) 


(WHz-') 





(10'" y) 

(mJy) 

(mJy) 

(10* K) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

J0953+2836 

70.2 

12.0 

18.8 

158 

-0.50 

25.9 

1.9 

1.94 

0.082 

0.028 

14-22 

4.35 

10.0 

1.2 

J1435+3131 

30.6 

7.8 

11.7 

324 

-0.72 

25.4 

2.2 

1.90 

0.098 

0.15 

1.3-2.0 

5.1 

5.5 

4.8 

J1722+5654 

22.8 

7.9 

12.2 

688 

-0.64 

25.3 

13 

1.52 

0.31 

0.20 

0.6-1.0 

5.4 

10.2 

9.8 


Note. — Columns are as follows: (1) lAU name (J2000); (2) projected full span of extended radio emission; (3)-(4) flux density of the core component and the 
total integrated flux density from the 9 GHz VLA data; (5) peak intensity of the 9 GHz emission in the approaching lobe; (6) non-simultaneous 1.4-9 GHz spectral 
index; (7) extended isotropic luminosity, extrapolated from 9 GHz to 1.4 GHz and A:-corrected assuming (Text = (8) ratio of peak intensity of approaching and 

receding lobes; (9) ratio of optical to radio core luminosity, /:-corrected to the emission frame assuming (Tcore = ‘^'ext = and (Topticai = -0.5; (10) projected jet 
speed computed from Rhs? (11) projected jet speed computed from jet/counter-jet length ratio; (12) kinematic age of source implied by (/?cos0)x for 6 = 10-15°; 
(13) 7.6 GHz VLBA core flux density; (14) total 7.6 GHz VLBA flux density; (15) emission frame VLBA core brightness temperature. 



A R.A. (arcsec) 


Fig. 1.— Naturally weighted maps of J0953+2836. Contours increase by 
factors of ^|2. At this redshift, I” = 6.97 kpc. (a) VLA map (9 GHz). To 
emphasize low-surface brightness extended features, a Gaussian taper with 
half-maximum at 1.5 x 10^ wavelengths has been applied. The first contour 
is at 50 pJybeam^^, about three times the map rms of 17 p,Jybeam“^, and 
the map peaks at 14.6 mJybeam“^. The elliptical Gaussian restoring beam is 
shown (FWHM 0.91 1" x 0.786^0- (b) VLBA map (7.6 GHz) of the central 
bright feature. The first contour is at 0.3 mJybeam“^, about three times the 
map rms of 0.1 mJybeam“^, and the map peaks at 5.7 mJybeam“^. The 
elliptical Gaussian restoring beam is shown (FWHM 1.82 mas x 1.39 mas) 

brightened hotspots near the outer edges of the kiloparsec- 
scale lobes evident at least on the approaching side. All three 
sources are moderately core dominated at these scales, with 
about 2/3 of the total emission in the kiloparsec core. The 



2 1 0 - 1-2 
A R.A. (arcsec) 


Fig. 2.— Naturally weighted maps of J1435+3131. Contours increase by 
factors o^^|2. At this redshift, I” = 6.12 kpc. (a) VLA map (9 GHz). 
The first contour is at 30 |Tjybeam“b about five times the map rms noise 
of 6 [iJybeam“b and the map peaks at 8.39 mJybeam“b The elliptical Gaus¬ 
sian restoring beam is shown (FWHM 0.236^' xO.228'0- (b) VLBA map 
(7.6 GHz) of the central bright feature. The first contour is at 0.3 mJy beam“b 
about three times the map rms noise of 0.1 mJybeam“b and the map peaks 
at 5.0 mJybeam“b The elliptical Gaussian restoring beam is shown (FWHM 
1.95 masx 1.24 mas). 

extended isotropic luminosities computed by subtract¬ 
ing the fitted core flux density from the total flux density and 
scaling to 1.4 GHz assuming a spectral index aext = -1, lie in 
the transition range between FR I and FR II galaxies, and fall 
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A R.A. (arcsec) 

Fig. 3.— Naturally weighted maps of J1722+5654. Contours increase by 
factors ofV2. At this redshift, 1" = 5.58 kpc. (a) VLA map (9 GHz). The 
first contour is at 21 pjybeam^’, about three times the map rms noise of 
7 pJybeam^’, and the map peaks at 8.14 mJybeam^*. The elliptical Gaus¬ 
sian restoring beam is shown (FWHM 0.341" X 0.221"). (b) VLBA map 
(7.6 GHz) of the central bright feature. The first contour is at 0.2 mJy beam^*, 
about three times the map rms noise of 0.07 mJybeam^*, and the map peaks 
at 4.8 mjybeam^'. The elliptical Gaussian restoring beam is shown (FWHM 
1.71 masx 1.21 mas). 

along the L-z trend reported in Kharb et al. (2010). 

On parsec scales, the sources are compatible with one¬ 
sided jet morphologies. The parsec-scale core region in 
J1435H-3131 is especially blazar-like, consisting only of a 
compact core. The other two sources are compatible with an 
alternative interpretation as a faint or undetected core with 
two-sided emission. 


4.2. Viewing Angles 

Two lines of reasoning suggest that the jet axes in these 
sources are at relatively small angles to the line of sight. First, 
the extended parsec-scale jets, where detected, are aligned 
with the kiloparsec-scale jets (this would not be true for 
J1722-b5654 if it is misaligned by about 180°). Misalign¬ 
ments due to projection effects are expected to be more com¬ 
mon for small viewing angles, so the alignment we find sug¬ 
gests the viewing angles are not especially small, although our 
statistics are limited (Pearson & Readhead 1988). 

Second, we reach a similar conclusion from R\, the ratio 
of the observed radio core luminosity to the optical luminos¬ 
ity. The latter acts as a proxy for the unbeamed, intrinsic 
emission. Motivated by Wills & Brotherton (1995), Kharb 
et al. (2010) found Ry to be a good indicator of viewing an¬ 
gle, where logRy = logLcore - logMabs - 13.69 with Lcore the 
1.4 GHz core luminosity and Mabs the k-corrected optical V- 
band magnitude. Assuming a flat core spectral index, we find 
that the logRy values for these NLSls (1.5-1.9) are at the low 
end of the distribution for jet-selected radio-loud AGN (e.g., 
the MOJAVE sample; Lister et al. 2013; Kharb et al. 2010). 
This suggests the viewing angles of these NLSls are similar 
to the highest angles in the MOJAVE sources, roughly 10- 
15° (Savolainen et al. 2010). 

4.3. Jet Speeds 


Assuming the structure and evolution of these sources is in¬ 
trinsically symmetric, we can estimate the projected jet speed 
in two ways. This common assumption is supported in these 
particular cases by the alignment of the approaching parsec- 
scale jet and the longer kiloparsec-scale lobe (Arshakian & 
Longair 2000). Both methods have a degeneracy between the 
speed and the viewing angle 0, but assuming 6> 10°, as seems 
to be the case, /3cos 0 ~ /3 to within the likely accuracy of our 
measurements. Eirst, we consider the flux ratio between the 
hotspots in the approaching and receding jets. If absorption is 
negligible along the additional optical path to the counter-jet, 
the expected ratio is 


/ 1 -H(j8cOS0)hs \^ " 
\1 -(y6COS0)hs/ 


( 1 ) 


where (J3cos 0)hs is the projected speed at which the lobes ad¬ 
vance at angle 0 to the line of sight. Eor this calculation, we 
use the ratio of the peak intensity (in pjy beam"*) between the 
jet and counter-jet lobes, (Rhs = We obtain sim¬ 

ilar results if we instead use the ratio between the flux den¬ 
sities obtained by fitting a Gaussian component model. Eor 
J0953H-2836 and J1722-1-5654, if we interpret the parsec-scale 
jets as two-sided, the flux density ratios of 1.2 and 2.3 imply 
jet speeds jScos0=O.O4 and 0.18, comparable to those found 
at kiloparsec scales. 

We can also examine the asymmetry in length between the 
jet and counter-jet. If the structure is intrinsically symmetric, 
differential light travel time will cause the approaching jet to 
appear longer than the receding jet (Banhatti 1980). In this 
case, the ratio x={q- rcj)/(rj + rq) = (J5cos0)x where rj and 
Tq are the lengths of the jet and counter-jet, and (J3cos0)x is 
the jet speed estimate. This method gives comparable results 
to the beaming ratio method, so overall we find/lcos 0 ~ 0.03- 
0.3. These correspond to Doppler factors in the range 5=1- 
1.5, where <5 = (1 -y6^)*^^/(l -/3cos0). 

Doppler beaming will increase the apparent brightness tem¬ 
perature of an emitting region relative to its intrinsic bright¬ 
ness temperature, with the observed Tq = ^TB^int. Us¬ 

ing our VLBA data, we computed Tq in the emission frame 
from an elliptical Gaussian model as Tb - (1.22 x 10*^ K)x 
5(1 + z)/(amajaminV^) where v is the observing frequency in 
GHz and 5, Omaj, and Omin are the flux density in janskys and 
the major and minor axis EWHMs (in mas) of the core compo¬ 
nent. The brightness temperatures we find are all substantially 
below the equipartition limit of about 10** K, compatible with 
a low Doppler factor (Readhead 1994). In each source, the 
presumed core has the highest brightness temperature of all 
fitted components by at least a factor of eight. 

The kinematic ages of these sources, computed assuming 
the approaching hotspots advanced outward steadily at their 
current speeds (/? ~ 0.03 - 0.3) at a viewing angle 0 ~ 10- 
15°, are found to be a few times 10^-10^ y. These are in 
the range found for other extended radio-loud NLSls, con¬ 
sistent with NLSls being much younger than blazars or radio 
galaxies (Doi et al. 2012). 


5. CONCLUSIONS 

Extended emission on kiloparsec scales was detected in all 
three of our targets. This is perhaps surprising in light of the 
small number of spatially extended NLSls previously known. 
However, few have been observed with the requisite angular 
resolution and sensitivity to detect the extensions we found 
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here. Because our targets were arbitrarily selected from a par¬ 
ent sample of bright, radio-loud NLSls, there was no reason 
all three were especially likely to show extended emission. 
Thus, our findings suggest that other NLS Is in the parent sam¬ 
ple are also likely to be extended. Observations of a larger 
sample are needed to verify this and, if true, to investigate 
whether this is generally true of all radio-loud NLSls or is a 
peculiar property of the parent sample. The morphologies of 
these sources on both kiloparsec and parsec scales are con¬ 
sistent with the presence of a mildly relativistic jet at a mod¬ 
erately small (10-15°) angle to the line of sight. The three 
sources we observed and most previously known extended 
NLSls where resolution permitted classification show FR II- 
like edge-brightened lobes. In contrast, Doi et al. (2014) re¬ 
cently reported edge-darkened FR I-like extended radio emis¬ 
sion in the NLSl Mrk 1239. The latter is radio quiet, while 
the others are radio loud, suggesting that radio loudness may 
be connected to intrinsic source properties, such as jet power, 
rather than orientation effects. 
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